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Abstract 

We report the measurement of longitudinal electron diffusion coefficients in liquid argon for electric fields 
between 100 and 2000 V/cm with a gold photocathode as a bright electron source. The measurement 
principle, apparatus, and data analysis are described. In the region between 100 and 350 V/cm, our results 
show a discrepancy with the previous measurement [1]. In the region between 350 and 2000 V/cm, our 
results represent the world’s best measurement. Over the entire measured electric field range, our results 
are systematically higher than the calculation of Atrazhev-Timoshkin [5]. The quantum efficiency of the 
gold photocathode, the drift velocity and longitudinal diffusion coefficients in gas argon are also presented. 

Keywords: Liquid Argon Time Projection Chamber, Electron Diffusion, Electron Drift Velocity, 
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1. Introduction 

The development of the liquid argon (LAr) ionization chamber, pioneered by Willis and Radeka [3] was an 
important advance in high energy calorimetry, and was quickly followed by the development of the LAr time 
projection chamber (TPC) [1] as a fine-grained tracking calorimeter for high energy physics experiments. 
LArTPCs are now the preferred technology for many accelerator neutrino and dark matter experiments. 
At present, two LArTPCs have been constructed and operated for neutrino physics measurements: the 
ICARUS [5] detector in Italy and the ArgoNeut [6] detector in the US. Meanwhile, two other LArTPCs 
have been constructed for dark matter searches: Darkside [7] and WArP [5] detectors. In the near term the 
MicroBooNE experiment, with a 170 ton LArTPC, has begun operation in the US In the future, a set 
of LArTPCs will be installed at Sanford Underground Research Facility (SURF) for the Deep Underground 
Neutrino Experiment (DUNE) to search for CP violation in the lepton sector and to determine the neutrino 
mass hierarchy [10] . For the near-term neutrino program m, a three-LArTPC configuration [12] will be 
implemented at Fermilab to search for a light sterile neutrino and to precisely measure neutrino-argon 
interaction cross sections. 

LArTPCs are attractive detectors for neutrino experiments. As the most abundant noble gas (^ 1.3% by 
weight) in the atmosphere, argon is commercially available in large quantities. The low cost and relative high 
density (~ 1.4 g/ml at 87 K) make LAr an ideal material for the massive TPCs needed for neutrino-induced 
rare processes. Their resulting charged particles transverse through the LAr and produce ionization electrons 
and scintillation light. Electrons in the ionization track will then drift at constant velocity along the lines 
of an applied electric field. The coordinate information perpendicular to the electron drift direction can 
be determined with a high-resolution, two-dimensional charge detector (e.g. wire planes). The scintillation 
light provides a fast indication of initial activity. By measuring the time delay to the subsequent signal from 
the drifting charge, it is possible to determine the distance over which the electrons drifted. The spatial 
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resolution of the interaction point and the detailed topology of the subsequent particle trajectories can reach 
the millimeter level. The reconstructed event topology can be used for the particle identification and, along 
with measuring the total amount of drifted charge, it can be used to determine the energy deposited in 
the detector. The charge collected per wire per unit time is closely related to the energy deposition per 
unit distance (dE/dx), which can also be used for the particle identification and the energy reconstruction. 
Therefore, the excellent signal efficiency and background rejection that result from these characteristics 
make LArTPCs ideal for neutrino experiments. 

In order to fully optimize the extraction of the intrinsic physics information from the recorded charge 
signal and to properly simulate the performance of LArTPCs, knowledge of the transport properties of 
electrons in LAr is essential. In particular, the diffusion of electrons drifting in the electric field from the 
point of ionization to the anode read-out plane is an important quantity contributing to the ultimate spatial 
resolution of the future long-drift-distance (up to 20 meters) detectors. The diffusion of electrons in strong 
electric fields is generally not isotropic. Therefore, longitudinal and transverse diffusion require separate 
measurements. For most substances the diffusion in the direction of the drift field (longitudinal diffusion) is 
smaller than the diffusion in the direction transverse to the field (transverse diffusion). Measurements have 
been reported previously of transverse diffusion at electric fields above 1500 V/cm [T311I1] and for longitudinal 
diffusion between 100 and 350 V/cm [J. In this paper we report a complete set of measurements of the 
longitudinal diffusion coefficient for electrons drifting in LAr at fields from 100 and 2000 V/cm. In Sec.[^ we 
review the basic formalism of electron diffusion in LAr. In Sec. we describe our experimental apparatus. 
In Sec. 1^ we describe the procedure of data taking and the analysis of the raw waveforms. In Sec. we 
discuss the systematic uncertainties in our measurement. In Sec. we report the results of electron drift 
velocity and diffusion in liquid argon as well as gas argon. A summary is presented in Sec. 


2. Electron Diffusion in Liquid Argon 


When a macroscopic swarm of electrons moves through a medium under the influence of an electric field, 
three processes are necessary to adequately describe its time development. They are 1) the drift velocity of 
the swarm centroid, 2) the diffusional growth of the volume of the swarm, and 3) the loss or gain of electrons 
in the swarm due to attachment to atoms or molecules in the medium or to ionization of the medium. In 
liquid argon at the fields considered here, ionization does not occur. The differential equation describing the 
time evolution of the electron density in the swarm is expressed by (Pick’s equation [T5l [T6] l: 
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where n = n(x,y,z,t) is the electron charge density distribution at position {x,y,z) and time t. The drift 
occurs with velocity v in the z (longitudinal) direction as shown in Fig. r; is the drift velocity. and 
Dt are the longitudinal and transverse diffusion coefficients, respectively. The attachment coefficient minus 
the ionization coefficient is represented by A. Since ionization does not contribute in our case, the inverse of 
A is equal to the mean free path. The solution to this equation for an initial point source of charge at the 
origin with a constant field in the 2 ; direction is described by a distribution function 
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where -\- y^ is the transverse coordinate. This distribution is a Gaussian function of 2 : and p at an 

instant in time, but is not a Gaussian in time at a fixed point in space. The electron current measured on 
a plane perpendicular to the drift direction at a distance d from the point source (i.e. the anode) is 
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This function approaches a Gaussian distribution for large d ■ v and small D ^, which is a reasonable approx¬ 
imation for our study. For this signal, the time at the peak, tp, is given by: 
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where d is the drift distance. The measured drift velocity vm is defined as the drift distance d divided by 
the time at the peak tp. The longitudinal diffusion width in time ctl, defined as the standard deviation in 
time of this charge distribution, can be similarly computed in terms of the true drift velocity v, the diffusion 
coefficient and the attachment coefficient A. The resulting expression is 
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with Vc = \/v{v + ADl\). For the longitudinal standard deviation we use the term diffusion time. 

The expression for the measured velocity and the diffusion time can be expanded in a power series in 
Dl to display the dependence on the parameters. The results of expanding both vm and ctl in powers of 
Dl are 
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Provided that the dimensionless variables a = DL/{dv) <^1 and /3 = adX <Cl, the values of vm and 
are well approximated by the first (zero order) terms in the expansions. Note that for attachment coefficients 
as large as A = 1/d, the condition that /3 <C 1 is still satisfied. The full expression for vm can be inverted to 
give V in terms of vm, Dl and A, but the expression for cr|^ cannot be inverted. However, by substituting 
the expression for v into the second series above, and applying Lagrange reversion m to the result, we can 
obtain the series expansions for the true quantities in terms of the measured ones. To simplify the result, 
we define the measured diffusion coefficient as 
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In terms of the directly measured quantities with vm = d/tp, it is 
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With this definition, the two inverted series can be expressed in terms of the small dimensionless variables 
ctM = Dl,m/{ dvM) and Pm = aMdX as: 
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These expressions are used to obtain the true drift velocity and diffusion coefficient from the measured 
quantities. 

We will show in Sec. that retaining only the first term in both these expansions results in a negligible 
bias for the measurements reported here. Note that the attenuation coefficient does not appear in the 
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expression for the diffusion constant: electron attachment has no effect on the measurement of the diffusion 
coefficient. Even for the drift velocity, the effect of electron attachment is still negligible: an attenuation of 
the electron swarm as large as 60% (Ad=l) will only triple the already small error that occurs when a <C 1 
(compare the second and fifth terms in the series for v in Eq. (§)■ The diffusion time, in terms of the true 
quantities, can be expressed to the same approximation as 

+ (^ - 3dA)a + ^(-313 - 680dA + + 0{a^)^ . (9) 

Again, for most purposes, only the first term needs to be retained. 

In these experiments we have recorded data directly from the output of the charge-sensitive pre-amplifier, 
so the measured signal is the total charge collected on the anode of the drift cell as a function of time, which 
is the integral of Eq. The signal shape is then described approximately by an error function, and the 
drift time and diffusion time are determined from the time that the signal reaches 50% of the maximum 
and half of the rise time between the 15.9% and 84.1% points of the signal, respectively. The drift velocity 
and the diffusion coefficients are functions of the applied electric held. To the extent that the higher order 
terms in Eq. are negligible, the longitudinal diffusion time is proportional to the square root of the drift 
distance at a hxed electric held. Fig. [T] illustrates the diffusion process. 




Figure 1: Illustration of the diffusion process starting from a point source to the detection plane indicated by dashed line. The 
diffusion in the longitudinal direction (along the drift direction) is generally expected to be smaller than that of the transverse 
direction. The standard deviation of the electron swarm (indicated by the shaded ellipses in the figure at three drift times) is 
expected to follow Eq. §. 

As the held approaches zero, the electrons gain so little energy from the held between the elastic atomic 
collisions that they come to thermal equilibrium with the medium. In this case, the dihusion coefficients are 
given by the Einstein-Smoluchowski relation [Hill]: 

hT 

D = —fi, (10) 

e 

where kT = 0.0075 eV for argon at the normal boiling point (87.3 K), e is the electron charge, and fjL 
is the electron mobility (the electron drift velocity v per unit electric held E). At low helds, when the 
electrons are in equilibrium with the medium, the mobility and the electron temperature become constants. 
As the held is increased and the electrons are no longer in thermal equilibrium with the medium, the 
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Einstein-Smoluchowski relation is taken to define the electron temperature. The electrons acquire energy 
from the field so that their temperature exceeds the thermal limit. The mobility decreases and the diffusion 
coefhcients increase. This process has been quantified by generalized Einstein relations which obtain a 
quantitative relation between diffusion and mobility for electrons in strong electric fields using the energy 
dependence of the electron-atom collision cross section. The earliest version of this class of relations for 
electrons, as proposed by Wannier nnun], has the form: 


n 

Dt = —/i, 
e 

kT, du . 


( 11 ) 


where kT is the electron temperature. The ratio of the longitudinal to the transverse diffusion coefficient is 
then expressed as: 


Dj^ dE 


( 12 ) 


At low electric fields, the mobility approaches a constant, and we have Dl/Dt « 1. As electric field 
increases, if the mobility decreases, then < Dx, and vice versa. 


3. Experimental Apparatus 

Historically, there are three different techniques for the measurement of longitudinal diffusion: un- 
gridded gridded [Mj, and shuttered |2S] drift cells. The electron source has typically been radioactive 

or a photocathode driven by a ultraviolet (UV) lamp. In this measurement, we chose to use a gridded drift 
cell because of its simplicity and similarity to the implementation of typical LArTPCs. For electron source 
generation, we employed an UV pulsed laser driving a semitransparent metal photocathode. The short pulse 
width and trigger stability of the laser allow timing information to be precisely synchronized. 



Figure 2: (color online) Illustration of setup for the electron drift velocity and diffusion measurement. See text for further 
discussion. 
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Fig. shows a schematic view of the experimental apparatus. The photocathode is composed of a 22 nm 
thick Au film 50% UV transmission) thermally evaporated onto a 1 mm thick, 10 mm diameter sapphire 
disk. The sapphire disk is mounted on a GIO dielectric ring at the top of the drift stack. The laser source 
is a frequency quadrupled, 4.66 eV, 60 picosecond pulse width Nd:YAG laser operating at 10 Hz repetition 
rate. The few /rJ of pulsed UV photons are focused onto a 1-meter long 600 p,m diameter UV fused silica 
fiber with an numerical aperture (NA) of 0.22. The fiber is then coupled to a 600 /rm diameter vacuum fiber 
feed-through followed by another 0.35-meter long 600 /im diameter in-vacuum UV fiber. The latter fiber 
ends directly on the back surface of the sapphire substrate. Photoelectrons are extracted into the LAr from 
the surface of the Au photocathode which is back-illuminated through this optical fiber delivery system with 
a maximum full beam size of ~ 1 mm. The laser energy density at the photocathode is at least an order of 
magnitude lower than the measured ^ 6 mJ/cm^ threshold which leads to damages of gold films. 

Negative high voltage is connected to a 200 ^m thick phosphor bronze disk gasket that makes electrical 
contact to the photocathode. Ceramic screws are used to secure the gasket onto the GIO disk structure. 
Field shaping copper electrode rings with 2 mm thickness, 21.3 mm outside diameter, and a 15.7 mm open 
aperture are spaced 5 mm apart, on center, dictated by four GIO precision slotted rod structures. A series of 
1 GH cryogenic compatible resistors make electrical contact to each ring and to the photocathode. A 300 /im 
thick Kapton ring is sandwiched between a -1-300 V biased anode disk below a grid mesh which is connected 
to the ground and to the last 1 GH resistor. The grid mesh has a wire width of 35 /tm with 350 /tm spacing 
and is mounted below a 250 /im thick phosphor bronze ring. The impact of thermal contraction on the 
distances at cryogenic temperature is negligible due to low thermal expansion coefficients of the materials. 
Given this setup, various drift distances from 5 mm to 60 mm can be arranged by adding or removing field 
shaping rings and moving the cathode assembly. The field between the grid mesh and the anode is called 
the anode collection field and kept at constant high value independent of the field in the drift region (drift 
field). 

A thermocouple is placed near the photocathode to monitor the temperature. The drift stack is housed 
in a 0.5 liter cylindrical vacuum cell capable of being pressurized up to 1.66 bar limited by a pressure relief 
valve. As shown in Fig.[^ electron charge arriving at the anode plane is fed to a charge sensitive pre-amplifier 
(10535) developed by Brookhaven National Laboratory which resides outside the cryogenic cell in ambient 
conditions. The pre-amplifier charge signal and its corresponding shaped and amplified pulses are recorded 
on an oscilloscope triggered by a fast photodiode that views a portion of the laser beam. The number 
of photoelectrons per pulse drifting through the stack to the anode ranges from 10^ to 10^ at bias fields 
increasing from 0.01 to ~ 3 kV/cm. On the low end, this measurement is limited by amplifier noise which is 
~ 1350 rms electrons. With an additional 0.1 /rs shaping amplifier, the noise charge can be further reduced 
to ~ 230 rms electrons. In addition, the pre-amplifier responds to noise induced by the laser Q-switch 
particularly near the t = 0 trigger time. Therefore, in our procedure, 64 to 512 signal traces are averaged 
and subtracted from the corresponding signal trace acquired when the anode extraction field is off. This 
data acquisition process takes about one minute per measurement at a given drift field. 

The entire system is baked to ^ 90°G and pumped to ^ 10“® Torr prior to measurements. Great 
attention is taken to achieve the necessary high purity of argon gas and to avoid discharge at the high bias 
fields [26]. We start with high purity commercial (99.9999 vol %) Ar gas, which is then passed through a dry 
ice cold trap to reduce water vapor. The Ar gas is further purified by an Oxisorb filter to a ppb impurity level, 
and then transported through stainless steel pipes and valves. The argon in the cell is not further purified by 
recirculation. Finally, Ar is liquefied by external cooling with a mixture of liquid nitrogen and dry ice. A CCD 
camera is mounted on a top viewport to observe the LAr liquefaction process and also serves as a liquid level 
monitor. To measure the work function and the quantum efficiency of the photocathode, a fiber coupled 
continuous wave UV white light source (Energetiq LDLS) is wavelength selected by a monochromator, 
fiber-coupled to the in-vacuum fiber feedthrough of the LAr drift cell, and the photocurrrent leaving the 
photocathode is measured by a Keithley picoammeter connected to the photocathode. 
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Figure 3: (color online) Summary of all data taken in liquid argon (LAr) and gas argon (GAr) in terms of the drift distance 
and the electric field. Data above 0.1 kV/cm was used for diffusion measurements. The boundary at high fields for each drift 
distance is limited by high voltage breakdown. 


4. Measurement and Raw Waveform Analysis 


Electron drift velocity and longitudinal diffusion are determined with data taken at various drift distances 
and electric fields. This is crucial to fully categorize systematic uncertainties. Changing the drift distances 
requires breaking the vacuum and rearranging the drift stack. Fig. shows the measurement parameter 
space in terms of the drift distance vs. electric field. At each drift distance, the maximum electric field that 
can be applied while avoiding high voltage breakdown is ~ 3 kV. It is expected that this limit may be raised 
with an apparatus design that reduces sharp points causing field emission. 

For each drift distance, the data taking procedure is 

1. measurement of electron charge versus drift field at 10“® Torr vacuum at room temperature, 

2. injection of ultrahigh purity Ar gas at up to 1.5 bar pressure into the test cell, following by charge 
measurement in gas Ar at room temperature, 

3. conductive cooling of the drift cell to the LAr temperature while passively maintaining the temperature 
at 87 ± 1 K and pressure of the drift cell to 1.2 ±0.1 bar of gas Ar, 

4. charge measurements at various drift fields in LAr. 


The cooling process in step 3 generally takes ~ 1 hour; the measurements in step 4 typically take several 
hours. It is important to note that LAr drift data are taken only when the thermocouple reading is 87 ± 1 
K and the CCD camera indicates no obvious turbulence in the liquid argon. 

The pressure measurements imply a temperature ranging from 88.2 K to 89.8 K at the surface of the 
LAr, in coexistence with the gas. This is distinctly higher than the temperature of 87 ± 1 K measured at 
the cathode. We attribute this difference to a presumed temperature gradient established in the LAr by 
cooling the Ar container from the bottom with the top near room temperature. Since no solid formation 
was observed during the measurements, we assume a minimum temperature of 84 K at the bottom of the 
cell. For these reasons we assign a mean temperature for the measurements in LAr of 87^3 K. The effect of 


the temperature and pressure uncertainties is discussed in Sec. 6.2 and Sec. 6.3 


Fig. [4] shows a representative set of pre-amplifier output data taken in LAr at 20 mm drift distance. 
Time zero is defined by the laser pulse trigger. The rapid rise of the signal corresponds to the time when 
the electron swarm generated at photocathode has passed through the grid mesh. The position and width 
of the rise correspond to the drift time and diffusion time, respectively. The slow exponential signal decay 
is due to the discharge by the feedback resistor on the pre-amplifier. The slow rise prior to the onset of 
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20 mm Drift Distance in LAr 



Figure 4: (color online) Representative raw signals in LAr at 20 mm drift distance are shown for various electric fields. The 
rising edges of the signals migrate toward the time zero with the increase of the electric field. 


the main signal is due to some non-negligible residual charge leaving the anode or arriving at the cathode 
and does not alter the drift electron signal shape. The reduction of signal height as the electric field is 
decreased due to the decrease in quantum efficiency at low electric field and to the finite electron lifetime in 
LAr. The effect of the electron space charge is negligible because of the low charge density produced at the 
photocathode (see Ref. [SI), which will be further discussed in Sec. 

The raw signal is fitted with a functional form of 

a, td, a) =a- e • ^1 -b Erf ( ^ 

which is a convolution of a Gaussian function approximating the electron swarm and a step response function 
with damping representing the response of the pre-amplifier. Here, a is the normalization factor that is 
proportional to the number of electrons in the swarm, a is the standard deviation of the Gaussian function 
which describes the longitudinal diffusion, rn is the decay constant of damping, and t is the drift time. To 
accommodate the slow rise before the main signal, an additional parameter for the baseline shift is added. 
A typical fit is shown in Fig. The uncertainty of each data point is the size of electronic noise, which 
is estimated by the standard deviation of the baseline. This type of uncertainty is treated as uncorrelated 
during the fit, and statistical uncertainties of the fitted parameters are adjusted accordingly to make the 
reduced chi-square x^/NDF unity. Systematic uncertainty due to the finite fitting range is estimated by 
changing the fitting range by 20%. 

5. Systematic Uncertainties 

The systematic errors in measuring drift velocities and diffusion coefficients caused by neglecting the 
higher order terms in the expansions of Eq. ([^ are negligible for our measurements. We have measured the 
electron lifetime in our apparatus, and found that it ranges from 66 (A = 0.3 cm“^) to 530 /is (A = 0.009 
cm“^). For these conditions, the expansion variables DL/{dv) and D^X/v are at their maximum values at 
lowest fields and shortest drift distances. For 5 mm drift distance and 0.1 kV/cm field, they have values of 
2.2 xlO”"^ and 3.4 xl0“^ respectively. The largest bias in neglecting all but the first term in the expansion 
of Eq. (|^ for this field and drift distance are 0.03% for drift velocity and 0.3% for diffusion time. For fields 
at which the electrons are in thermal equilibrium with the LAr ( 0.1 kV/cm), the value of a can be 

computed using the Einstein-Smoluchowski relation to be a = cl!, where AU is the potential drop over 
the drift distance and cl is the electron energy in volts. For higher fields and longer drift distances the 
errors are smaller. 
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Example Fit: 20 mm, E = 0.1 kV/cm 



Time (|j,s) 

Figure 5: An example fit of Eq. in is shown. The uncertainty of each data point represents the size of estimated electronic 
noise and is treated to be uncorrelated in the fit. The uncertainty of the data are enlarged to make the reduced chi-square 
/NDF to unity, which leads to increased uncertainties of fitted parameters. 


The second systematic uncertainty is from the impulse response function of the pre-amplifier. The 
integration time of the pre-amplifier is tj. In this case, the full form of signal response function is 


f{t,cr,TD,a,Ti)=a-e 1 + Erf 
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where the first term has the same functional form as that in Eq. (13). Due to the smallness of the t/, the 


above equation is difficult to evaluate numerically without a computational demanding integration. Since 


the full functional form is not suitable in the large-scale fitting of the data, Eq. (131 is used instead to extract 


the drift and diffusion time. In this case, t/ contributes to both the extracted drift and diffusion time (i.e. 
the rise time) of the electron signal by adding additional time, which needs to be corrected. This effect was 
directly calibrated by fitting the electron signal measured in vacuum. The contribution to the drift time 
and diffusion time due to t/ is determined to be 31 ± 2 ns and 33 ± 5 ns, respectively. The systematic 
uncertainty due to different fitting functions (Eq. (13) vs. Eq. 0) is negligible. 

The largest systematic uncertainty is caused by the electron drifting from the grid mesh to the anode. 
The grid mesh is installed at 300 /im from the anode. The purpose of the grid mesh is to screen the anode 
from the slow signal induced by the electrons in the drift region. The voltage at the anode is maintained at 
300 V and the grid is at ground, creating a 10 kV/cm electric field in the collection region, which is much 
higher than applied electric field in the drift region (0.01 - 3 kV/cm). A perfect grid would separate the 
drift and collection regions with a constant field on both sides. In reality, electric field “leaks” through the 
grid from the high-held region (i.e. the collection region) to the low-held region (i.e. the drift region). This 
effect leads to a bias in estimating the drift held. In order to evaluate this effect, the electric held in the 
drift cell is computed by MAXWELL 3D m using a geometric model of the apparatus. The scheme and 
resulting electron trajectories of the simulation are shown in Fig. Fig.j^shows the computed electric held 
for a set of points along the diagonal axis of the cell. The electric held is largely non-uniform within -200 to 
4-400 fim distance from the grid mesh. 

Another possible systematic uncertainty is the imperfect uniformity of the electric held in our drift 
stack. We used MAXWELL 3D to compute the electric helds for the full geometry including the drift 
stack, cryostat walls, and cables etc. Ray tracing calculations in these helds show a very small focusing of 
the electron swarm as it drifts to the grid. This causes the transverse width of the beam distribution to 
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Figure 6: (Left) The scheme of the electron signal simulation. The electric field is computed with the MAXWELL 3D in the 
entire geometry (i.e. including both the drift and collection regions) with proper boundary conditions. (Right) The simulated 
electron trajectories are plotted with different colors respect to the grid mesh. The width of a grid cell is 350 pm. The curved 
trajectories lead to a bias on the extraction of both drift time and diffusion time which are corrected in the analysis. See text 
for further discussions. 


decrease by about 0.08 mm, which is significant compared to the estimated transverse diffusion especially 
for short drift distances. However, the variation in the drift time caused by the variation in the path lengths 
is much smaller fractionally, systematically increasing the longitudinal diffusion time by a maximum of 0.31 
ns compared to the measured diffusion time of 42 ns for the highest fields reported. This is a maximum 
systematic increase of 0.7% in the longitudinal diffusion coefficient. This systematic, one sided error in the 
diffusion is typical of the calculations for all the drift lengths. Any other instrumental effects producing time 
invariant fields also can only increase the apparent measurement of diffusion. Among such effects that have 
not been modeled are: focusing through the grid mesh; non-co-planarity of the grid, anode, cathode and 
field rings; imperfect spacing and shape of the gradient rings; charging of support insulators; and imperfect 
matching of the gradient resistor values. Including the effects of these imperfections on the central field in 
the drift cell increases the assigned error for field imperfections to t.^%. This one-sided systematic error is 
still small compared to the statistical error in the analysis. 

The electron signals are simulated with a ray tracing technique that begins with randomly distributed 
electrons within the dimension of a single grid cell (350 /rm wide on each dimension). As shown in the right 
panel of Fig. 200 electrons with zero initial velocity are generated 5 mm above the grid plane. The electron 
trajectories are solved using the computed electric field and expected electron mobility from a global fit of 
the existing world data. The electron trajectories are focused into the center of the grid cell when electrons 
are close to the grid mesh, as shown in Fig. Therefore, three sources of systematic uncertainties arise due 
to the imperfect grid plane: 

• The drift distance is measured to the grid mesh. The electrons take some additional time to drift far 
enough toward the anode to cause the signal rise of 50% (the stop for the drift time measurement 
clock). In addition, the trajectories starting far from the center of the cell are longer than the central 
trajectory due to the focusing as shown in the right panel of Fig. This variation also increases the 
drift time. This is negligible for all but the shortest drift distances and highest drift fields as will 
be discussed in Section |6.2[ The tightly focused bundle of trajectories for low drift fields lead to a 
spreading of the drift times over the grid cell, which increases the signal rise time over that due to 
longitudinal diffusion only. This effect is expected to be largely independent of the drift field. The ray 
tracing simulation suggests that the magnitude of longitudinal diffusion is comparable to the electron 
drift time from the grid plane to the anode plane. Therefore, we use the expected drift time from grid 
to the anode plane to estimate the contribution of non-uniform electron trajectory to the diffusion 
time and assign a conservative 100% uncertainty in the analysis. This effect is negligible for GAr, but 
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Figure 7: (color online) MAXWELL 3D computed electric field close to the grid mesh are shown for a set of points along the 
diagonal axis of a grid cell. For the expected electric field in the collection and drift region of 10 and 0.1 kV/cm, the achieved 
fields are respectively lower and higher than the expected fields due to the field leakage through the grid. These curves are 
labeled with numeric indexes by group on the diagonal line. 
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Figure 8: (color online) The square-root of photocurrents normalized to the value at 5.16 eV from photocathode are shown 
with respect to the photon energy. Straight line fits on the data between 4.43 and 5.14 eV are performed to extract the work 
function at y=0 intercept. Representative error bars are shown in the fitted region. 



important for the LAr measurement. 

• When the drift held is small, the leakage held leads to a signihcant bias on the magnitude of the held 
in the drift region. As shown in Sec.|^ we limit our diffusion measurement to helds above 0.1 kV/cm. 
We further estimate the uncertainty due to this effect with an empirical formula (see Section [6(3| ) . 

Additional systematic uncertainties include the measurement of drift distance (measured to 0.1 mm) and 
the control of experimental conditions (pressure and temperature). For LAr the temperature is accurate to 1 
K. For GAr, it is controlled to the room temperature, and the pressure is controlled to 2 mbar. The resulting 
uncertainties will be discussed in Sec. |6.2| and Sec. |6.3| for electron mobility and diffusion measurements, 
respectively. 

6. Results 

6.1. Quantum Efficiency of Au Photocathode 

Several batches of Au hlms are measured and consistently have a work function of 4.17 eV, 4.10 eV, and 
3.89 eV in vacuum, in GAr, and in LAr, respectively, see Fig. The 0.28 eV lowering of work function in 
LAr indicates a negative electron affinity, which is similar to previous measurements with Zn |28j and Ba 
[35] photocathodes. Although the drift photoelectrons bear this excess kinetic energy on ejection from the 
cathode, they come to thermal equilibrium through collisions with argon atoms in a time less than 1 ns, and 
a distance less than a few microns [30] . Therefore the negative electron affinity has no impact on drift and 
diffusion measurement. 

The quantum efficiency (QE) is defined as the number of electrons leaving the photocathode per number 
of UV photons irradiating the back surface of the photocathode. QEs of ^ 10“®, ^ 2 x 10“®, and ~ 10“^ 
are obtained when the photocathode is operated in vacuum at room temperature, 1500 ± 2 mbar of high 
purity Ar gas at room temperature, and at LAr temperature, respectively. Despite the QE drop of nearly 2 
orders of magnitude in LAr, the photocathode provides 10® electrons with ^ 1/iJ of UV photons which are 
sufficient for the measurements described here. 

The factor of ^ 5 drop in the QE from vacuum to 1500 mbar of Ar gas is due to the increase of the 
potential barrier caused by the Ar atoms adsorbed on the photocathode, which can be recovered by pump 
down of the system and proper baking. The factor of ~ 50 drop in the QE from vacuum to LAr is mainly 
caused by electron back diffusion to the photocathode. Although the negative electron affinity of -0.28 eV 
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Figure 9: (color online) (Left panel) Charge density from the photocathode is plotted versus the electric field for vacuum, GAr 
and LAr. (Right panel) The square-root of the quantum efficiency is plotted against the square-root of the electric field. Linear 
fits of the data are shown as well. 


reduces the work function of the photocathode, the field of the image charge and the low drift velocity 
of electrons in LAr favor strong back diffusion so that most photoelectrons return to the photocathode. 
A random walk simulation of the diffusion process for an electron leaving the cathode with excess energy 
is presented in m- The results of this simulation suggest that when an electron starts on the gold-LAr 
interface with an electron energy ^ 0.5 eV and propagates in a held < 1 kV/cm, the survival rate of the 
photoelectrons is ~ 3 %. This result is consistent with the QE measurements, and also qualitatively agrees 
with other measurements on metal photocathodes operated in LAr [28] . 

Photoemission is a statistical process which can be described, to the hrst approximation, by the well- 
known three-step Spicer photoemission model [35]: (i) absorption of light to excite the electrons from an 
initial state to a hnal state, (ii) transport of excited electrons to the surface, and (hi) escape of the electrons 
by overcoming a potential barrier. The QE is thus given by the Fowler law: 

QE = k{hv — i/o + Cf \/P ■ EY: (15) 

where fc is a constant that depends on the material, (po is the work function of the metal, a = \/ e/47reo, e is 
the electron charge, Cq is the vacuum dielectric constant, E is the extraction electric field, hv is the photon 
energy, and /3 is the field enhancement factor. The left panel of Eig. plots the charge density leaving the 
photocathode versus the bias held, while the right panel of Fig. j^plots VQE against the The measured 
QE exhibits a linear dependence as predicted by Fowler’s law in vacuum and in low pressure of gaseous 
argon but has some deviations in liquid argon. A more quantitative evaluation of these deviations requires 
measurements outside the current scope. 

6.2. Drift Velocity in GAr and LAr 

The drift velocity and mobility are dehned as 


V = = fjL ■ E. (16) 

Here, d is the drift distance, t is the drift time, yi is the mobility, and E is the electric held. The drift distance 
d is measured to an uncertainty of 0.1 mm by the position of the anode and cathode on the drift stack. The 
drift time t is calculated from the measured drift time by correcting for the time shift due to the response of 
the pre-ampliher and for the time which electrons take to drift from the grid mesh to the anode plane. The 
pre-ampliher time offset is measured to be 31 ± 2 ns with respect to the laser trigger from the calibration 
data taken in vacuum. The drift time between the grid plane and the anode plane for GAr is calculated 
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Figure 10: (color online) Electron drift velocity (top panels) and electron mobility (bottom panels) are shown with respect to 
the electric field for various drift distance settings for GAr (left panels) and LAr (right panels). The nominal electric field E 
(top x-axis) and the reduced field E/N (bottom x-axis) which is the E normalized with the number density of argon with the 
unit of Td = 10“^^ V- cm^ are both presented for direct comparison. The red bands in the mobility plots (bottom panels) 
represent the global fit and their uncertainties. For LAr, the expected mobility at zero electric field (~ 518 cm^/V/s) is plotted. 
See text for further discussions. 
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to be 15 ± 5 ns based on drift velocity data from Refs. For LAr, this correction is calculated to 

65 ± 7 ns based on data from Refs. [29l |45H58]. The additional drift time due to grid-to-anode drift and 
curved trajectories is estimated with ray tracing simulation to be 15 ± 15 and 65 ± 65 ns for GAr and LAr, 
respectively. These corrections are small compared to the measured drift time ranging from a few /rs to a 
few hundred ^s. The uncertainties due to temperature and pressure variations are estimated to be 3% and 
5% assuming the general temperature dependence of the mobility [33]. For LAr, this uncertainty is 

confirmed by repeated sets of data taken over a period of several hours at 20 mm drift distance and 0.5 
kV/cm electric field. 

Beside these uncertainties, an additional systematic uncertainty is added due to the leakage field as 
discussed in Sec. The leakage field leads to a sizable difference in LAr electron mobility measured from 
different drift distances for electric field below 0.1 kV/cm. Due to the difficulty in precisely predicting the 
size of leakage field in all the configurations, no correction was made to our results. Instead, conservative 
uncertainties for GAr and LAr were estimated by quoting a 100% change in results when including empirical 
leakage field assumed to be uniform in the entire drift region: 


0.025 • 


5mm 


(17) 


in kV/cm. The format of the leakage field was then deduced, so that the zero-field electron mobility in LAr 
then satisfied the expected results of ~518 cm^/V/s. The impact of this effect is much smaller for GAr due 
to much weaker electric field dependence in our region of interest for the drift velocity. For the diffusion 
measurement, we limit the data to a region with expected electric field larger than 0.1 kV/cm, where the 
leakage field effect has a smaller effect. 


Fig. 10 shows electron drift velocity (top panels) and electron mobility (bottom panels) as a function of 


the drift field for various drift distances for GAr (left panels) and LAr (right panels). The nominal electric 
field E and reduced field E/N which is the E normalized with the number density of argon with the unit of 
Td = 10“^^ V- cm^ are both presented for direct comparison. This normalization is in particular importance 
for GAr comparison due to its wide range of density variation at different pressures and temperatures. The 
red curve and band in the bottom panels are the global fit of world data, in LAr data from Ref. |33J H5H55] , 
GAr data from Ref. . 

The relative uncertainties for the global fit are about 13% and 20% for LAr and GAr, respectively. 
For GAr the extracted electron velocities are largely consistent among different drift distances except for 
the 5 mm drift, which indicates additional experimental uncertainties unaccounted for in the analysis in 
this setting. As shown in the Sec. the drift time is a parameter used in the extraction of diffusion 
coefficient. A 7% additional uncertainty was added to the measured drift time in order to take into account 
this discrepancy. 

For LAr, the extracted electron velocities at various drift distances are consistent with each other. The 
measured velocities for 5 mm and 10 mm show a small discrepancy at high electric field, which again indicates 
additional experimental uncertainties unaccounted for in the analysis. A 5% additional uncertainty was 
added to the measured drift time in order to extract the diffusion coefficient conservatively. The comparison 
of the extracted LAr electron mobility data are in excellent agreement with the global fit for electric fields 
larger than 0.1 kV/cm. At electric fields smaller than 0.1 kV/cm, the leakage field effect leads to significant 
systematic uncertainties. The black line represents the expected electron mobility ~ 518 cm^/V/s at 87 K, 
which is used as a guideline to estimate the strength of leakage field. For the diffusion result discussed in 
the next section, we are restricted to drift fields larger than 0.1 kV/cm. 


6.3. Longitudinal Diffusion in GAr and LAr 

The longitudinal diffusion time defined in Eq. (|^ can be equivalently expressed as: 




2ei • d t 
E ' d' 


(18) 
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Figure 11: (color online) The extracted electron energy for the longitudinal diffusion is plotted versus the electric field for 
GAr and LAr. The nominal electric field E (top x-axis) and reduced field E/N (bottom x-axis) which is the E normalized with 
the number density of argon with the unit of Td = 10“^^ V- cm^ are both presented for a direct comparison. The red band 
in the GAr panel represents the global fit and its uncertainty. For LAr, our data are systematically higher than the prediction 
of Atrazhev-Timoshkin [2]- The red line represents the longitudinal electron energy and blue line represents the transverse 
electron energy. The magenta points are deduced from the measurements by IGARUS ^ at low electric field (see text for 
details). The green dots represent the transverse electron energy data reported by Shibamura m at high electric fields. The 
dashed cyan curve is the calculated longitudinal electron energy using the general Einstein relation (Eq. and Eq. ( |12[ |) 
from the interpolation (dashed green line) to the transverse electron energy reported by Shibamura m- The brown dot-dashed 
line represents the parameterization we introduced. Details are presented in the appendix. 
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Here, €l = cDl/^ is the effective energy of electron, Dl is the longitudinal diffusion coefficient, /x is the 
electron mobility, e is the electron charge, d is the drift distance, E is the electric field, and t is the drift time. 
In our experimental setup, we expect additional contributions from systematic errors due to the measured 
drift time. The contribution from the response of the pre-amplifier is determined with the calibration data 
in vacuum to be 33 ± 5 ns. As discussed in Sec. the contribution due to grid-to-anode drift and curved 
electron trajectories are estimated to be 15 ± 15 and 65 ± 65 ns for GAr and LAr, respectively. In this case. 


Eq. (18) is modified as: 


al = eUE) 


d-E 


\d), 


(19) 


with c{d) being the constant term that depends on the drift distance. The longitudinal diffusion time aL and 
drift time t depends on both the electric field E and the drift distance d. Therefore, the effective electron 
temperature eL{E), as a function of the electric field E can be extracted through a simultaneous fit of all 
the results of longitudinal diffusion and drift time by minimizing the following function: 


X' 


= E 


(ai - eL(Afc) • 


AaliSaL^ 


d^,Ei 
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{Sd,y 


iSE^)2 

E? 


4(sty 

*2 


E 


icjdj) - cqY 
(<5c)2 


( 20 ) 


where the notation SX represents the uncertainty of the physical quantity X. The notation i represents 
the ith data point with electric field E^ and drift distance dj. The constant term at each drift distance 
is constrained around the expected value with its uncertainty. Beside the uncertainties of measured drift 
time as discussed previously, the uncertainty of the electric held is estimated with the empirical formula in 
Eq. ( [T7| . Additional uncertainties of 1% and 3% are included for GAr and LAr, respectively, to take into 
account the variations in experimental conditions in terms of temperature and pressure. This estimation 
was obtained assuming a linear temperature dependence for the electron energy. Furthermore, due to the 
leakage held ehect, the ht was limited to data taken higher than 0.1 kV/cm drift held where the leakage 
held effect is properly estimated as suggested by the mobility data discussed in Sec. |6.2| 

For GAr, the reduced chisquare 


The extracted electron energy cl for GAr and LAr are shown in Fig. 11 


X^/NDE = 3.3/23 indicating conservative systematic uncertainties. The extracted electron temperatures 


are consistent with a global ht to the world data [37l 141114^ . The consistency between our data and the 
global ht again conhrms the validity of the analysis procedure. 

For LAr, the reduced chisquare x^/NDF = 114.6/43 indicating additional unknown systematic uncer¬ 
tainties. In order to minimize the bias to the results due to these uncertainties, an additional 8% uncorrelated 
uncertainty is added to each of the measured diffusion times. The new ht shown, in the right panel of Fig. [m 
gives a reduced chisquare /NDF = 45.8/43. The IGARUS group reports measurements for the square 
of the diffusion time vs. drift time at 92 K and four helds from 100 to 350 V/cm in Fig 15. of Ref. [T]. 
We have ht this data to obtain longitudinal diffusion coefficients at the four helds. They also report drift 
velocity measurements at 92 K for helds above 200 V/cm. In order to convert the diffusion constants to 
effective electron temperatures using the Einstein-Smoluchowski relation (Eq. @) , we have obtained the 
corresponding mobilities by htting their data for mobility at 92 K (contained in Fig. 13 of Ref. [1]) to the 
functional form we used for the global ht (Eq. (21) in the Appendix). These effective energies at 92 K are 


scaled to 87 K assuming a linear T dependence of the effective electron energy. This is certainly correct at 
low helds, where the electrons are in thermal equilibrium with the LAr, and is probably quite accurate for 
the entire range of the IGARUS data set, for which the mobility is still close to the zero held value. Our 
results are higher than the converted results from IGARUS. 

At high electric helds, there is no existing data on the longitudinal dihusion. On the other hand, 
transverse diffusion results exist for electric held > 1.7 kV/cm from Shibamura [13]. As discussed in Sec.[^ 
the longitudinal electron energy can be related to the transverse electron energy through the general Einstein 
relation. Using Eq. 0 and Eq. ( [T^ together with the electron mobility from the global ht (Fig. [Toj ), 
the transverse electron energy results from Ref. |13j (green points in Fig. 11 and green dashed line for 
interpolation) are converted into longitudinal electron energy (dashed cyan curve in Fig. 11). The deduced 
longitudinal electron energies are consistent with our data at high electric helds. 
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The formula of our fit to a large selection of the world’s data for electron mobility in LAr is presented 
in the Appendix. A similar parameterization is applied to the existing data on longitudinal electron energy 
as shown in Fig. El whose formula is also presented in the Appendix. 

Our results are also compared to a calculation of the transport properties of LAr by Atrazhev-Timoshkin [5], 
which gives the transverse and longitudinal electron energy as a function of field near the triple point (their 
Fig. 4). These calculations are shown in Fig. 11 Our data are systematically higher than this calculation. 


At zero electric field, the electron temperature is expected to be the same as that of the LAr. The fact that 
both ICARUS’ and our data are higher than this value (0.0075 eV at 87 K) at low electric field is puzzling. 
In the ICARUS paper [T], an additional correction to (and thus to the electron energy) caused by the 
space charge effect (Coulomb repulsion of the electron cloud) was quoted. This correction was estimated 
by using an approximate model described in Ref. m- The resulting electron energy from ICARUS after 
applying the correction is much lower than the LAr temperature. We performed the same calculation and 
found that the correction for our conditions is small due to the lower electron density in our measurements 
for both drift and collection region. The focusing of the electron trajectories through the grid mesh causes a 
negligible increase in any space charge broadening since the transit time through the grid mesh is a negligible 
fraction of the total drift time and the total charge is divided among many grids. We also note that Ref. m 
declined to apply space charge effect they computed to their measurements. We agree with Ref. m and 
do not include the space charge effect correction in our results. Nevertheless, we should comment that 
the space charge effect at low field with short drift distance and maximum space charge density could be 
more significant and may possibly explain the discrepancy between our effective energy to the theoretical 
calculation. This can best investigated by measurement with different space charge density, which can be 
achieved by changing the transverse source dimensions. 


7. Summary 

LArTPCs is a key detector technology for neutrino physics in searching for new CP violation and de¬ 
termining the mass hierarchy. Understanding the performance of LArTPCs requires a good knowledge of 
electron transport in LAr. In particular, diffusion of electrons in LArTPCs limits the spatial resolution of 
tracks with long drift distances. At 500 V/cm electric field, the longitudinal diffusion coefficient is Dl = 7.2 
cm^/s from measurements, and the transverse diffusion coefficient is Dt = 12.0 cm^/s by extrapolating 
existing data measured at high field [T3]. For a 3.6 meter drift, as proposed in the DUNE far detector, the 
diffusion limited resolution will be around 1.8 mm (1.1 ^s in time) longitudinal and 2.5 mm transverse to 
the electric field at 500 V/cm for a track produced near the cathode. This implies spatial resolution for 
the longest drift distance will not be improved with signal sampling time shorter than 1 /rs and wire pitch 
much less than 3 mm. For extremely long drift distances, such as the 20 meter drift distance proposed for 
the GLACIER detector, the diffusion limited resolution will be 4.2 mm (2.7 /iS in time) longitudinal and 
7.7 mm transverse to the electric field of 500 V/cm. These values become important for distinguishing very 
short tracks from point like depositions which must drift for long distances. An important example include 
distinguishing supernova neutrinos from the more nearly point-like tracks made by the a and /3 particles 
from the decay of radioactive contaminants in the LAr (such as ^®Ar). Since diffusion plays a non-negligible 
role in determining signal shape in large LArTPCs, a good knowledge of diffusion is important in simulating 
the device responses and in extracting track coordinates. 

In this paper, we present our measurement of the electron mobility and longitudinal diffusion coefficients 
in GAr and LAr. Our experimental setup uses a thin gold photocathode as a bright electron source. The 
work function and quantum efficiency of the gold photocathode are measured in vacuum, GAr, and LAr. 
Our result shows that the gold photocathode is a good candidate as a calibration source for future large 
LArTPCs. 

The measured electron mobility in GAr and LAr and longitudinal electron diffusion coefficient (expressed 
as the electron energy) in GAr, over a wide range of electric fields, show good agreement with fits of existing 
world data. This validates our experimental approach to systematically measure the longitudinal electron 
diffusion coefficients in LAr for electric fields between 100 and 2000 V/cm. The extracted longitudinal 
electron energy in LAr in the region between 100 and 350 V/cm shows a discrepancy with previous ICARUS 


18 



measurement. In the region of 350 to 2000 V/cm, our results represent the world’s best measurement. Over 
the entire electric held range, the measured longitudinal electron energy is found to be systematically higher 
than the prediction of Atrazhev-Timoshkin [2]. 

Our results are limited by the systematic uncertainties coming from i) leakage held ehects due to im¬ 
perfection of grid mesh and ii) unaccounted systematic uncertainties in the diffusion time likely due to the 
control of experimental condition and noise. A new experimental setup with improved design of the grid 
mesh and better control of experimental conditions has been constructed. Improved measurements on the 
electron temperature are expected and will be reported. 

It is also desired to extend the measurement to transverse diffusion. While the longitudinal diffusion 
expresses itself as a spread in drift time, the transverse diffusion expresses itself as a spread in space. There¬ 
fore, the anode in the current apparatus is not sufficient to perform measurements of transverse diffusion. 
A new anode with hne position resolution is being designed for future transverse diffusion measurements. 
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Appendix 

Global fit on electron mobility in LAr 

The electron mobility data set used in our global data ht are from Ref. [5SJ H5H55] . LAr temperature 
used in this global ht is 89 K. The data reported in the references are all scaled to this temperature with 
the temperature dependence of [^. The htting function is a rational polynomial expressed as: 

ao + aiE + 02^3/2 + (T\ 

^ l + {ai/ao)E + a4E^ + a5E3\ToJ ’ ^ ’ 

where E is the electric held in the unit of kV/cm and oq = 551.6 cm^/s is the electron mobility at zero held 
with temperature of Tq = 89 K, T is the LAr temperature. The htting parameters are given in Table 


Table 1: Function paraemeters in Eq. 113 

ao = 551.6 

ai = 7953.7 

02 = 4440.43 

03 = 4.29 

04 = 43.63 

05 = 0.2053 


Parameterization on ejfeetive longitudinal electron energy 

We also introduce a parameterization of the effective electron energy for the convenience of application. 
Both of our data and ICARUS’ data at low held are included. The parameterization is also in a form of 
rational polynomial 

^ bo+b^E + b^E^ / T \ 

^ l + {bi/bo)E + boE^\Tj ^ > 

where E is the electric held in the unit of kV/cm and bo = 0.0075 eV is the electron energy at Ti = 87 K 
under zero held and T is the LAr temperature. The parameterization can be applied to other temperatures 
with a linear temperature dependence of T. 
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Table 2: Function parameters in Eq. \22\ 

bo = 0.0075 

6 1 = 742.9 

62 = 3269.6 

63 = 31678.2 


The longitudinal diffusion coefficients Dl in the range of 0.1 to 1.5 kV/cm can thus be expressed as 

(23) 


defined by the Einstein relation in Eq. (10): 
ao 


D, = ^ = 

e 


. + aiE + 02^3/2 + a3E5/2 \ 


biE - 


b2E^ 


l + {ai/ao)E + a4E^ + a5Ey \1 + {bi/bo)E + boE^ J \T^ 
with the parameters given in Tables. and 
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